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ABSTRACT 

We discuss the influence of tfie cosmological background density field on the spherical 
infall model. The spherical infall model has been used in the Press-Schechter formalism 
to evaluate the number abundance of clusters of galaxies, as well as to determine the 
density parameter of the universe from the infalling flow. Therefore, the understand- 
ing of collapse dynamics play a key role for extracting the cosmological information. 
Here, we consider the modified version of the spherical infall model. We derive the 
mean field equations from the Newtonian fluid equations, in which the influence of 
cosmological background inhomogeneity is incorporated into the averaged quantities 
as the backreaction. By calculating the averaged quantities explicitly, we obtain the 
simple expressions and find that in case of the scale-free power spectrum, the density 
fluctuations with the negative spectral index make the infalling velocities slow. This 
suggests that we underestimate the density parameter fl when using the simple spheri- 
cal infall model. In cases with the index n > 0, the effect of background inhomogeneity 
could be negligible and the spherical infall model becomes the good approximation 
for the infalling fiows. We also present a realistic example with the cold dark matter 
power spectrum. There, the anisotropic random velocity leads to slowing down the 
mean infalling velocities. 
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1 INTRODUCTION 

In the standard scenario of the cosmic structure formation, 
the density peaks in the large-scale inhomogeneities grow 
due to the gravitational instability and they experience the 
gravitational collapse, which finally lead to the virialized 
objects such as the clusters of galaxies or the dark matter 
halos. 

Modeling such non-linear objects plays an important 
role for probing the large scale structure. After the seminal 
paper by Gunn & Gott (1972), the spherical infall model 
has been extensively used in the literature. The non- linear 
system mimicking the density peak is represented by the 
spherical overdensity with the radius R, which obeys the 
simple dynamics (Peebles 1980) 

(1) 

where M is the mass of the bound object. The density con- 
trast of this non-linear object is defined by 5 = {a/ Rf — 1, 
where a is the expansion factor of the universe. From (^, we 
can estimate the collapse time tc, which gives the dynamical 
time scale of the virialization. In terms of the linearly extrap- 
olated density contrast, we have the value 5c = S{tc) — 1.68 



in the Einstein-de Sitter universe (for the low density uni- 
verse or the flat universe with the non-vanishing cosmolog- 
ical constant, see e.g. Eke, Coles & Frenk 1996). This is 
frequently used in the Press-Schechter mass function to eval- 
uate the number abundance of cluster of galaxies (Press & 
Schechter 1974, Eke, Coles & Frenk 1996, Kitayama & Suto 
1997). 

Further, the spherical infall model can also be applied 
to estimate the density parameter of our universe, although 
the other efficient method such as POTENT has recently 
been available (e.g, Bertschinger & Dekel 1989, see also the 
review by Strauss & Willick 1995). The observation of the 
infalling velocity Vinf can be compared to the expansion rate 
R/R subtracting the uniform Hubble flows. Characterizing 
the infalling velocity Vinf as the function of mean overdensity 
5 within the sphere centered on the core of cluster, we can 
determine from the S-Vini relation (Davis & Huchra 1982, 
Yahil 1985, Dressier 1988). 

As we know, the spherical infall model is the ideal sys- 
tem. In real universe, the density peak cannot evolve lo- 
cally. The cosmological background density field affects the 
dynamical evolution and the collapsed time of the density 
peak. This might lead to the incorrect prediction for the 
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low mass objects using the Press-Schechter formalism (Jing 
1998, Sheth, Mo & Tormen 1999). As for the infall veloc- 
ity, the tidal effect distorts the flow field around the core of 
cluster, which becomes the source of the underestimation of 
the density parameter (Villumsen & Davis 1986). 

Therefore, the understanding of the influence of cos- 
mological background density field on the bound objects 
is crucial to clarify the cosmic structure formation. Sev- 
eral authors treat this issue and consider the modification 
of the spherical infall model taking into account the back- 
ground density inhomogeneity (Ifoffman 1986, Bernardeau 
1994, Eisenstein & Loeb 1995). The validity of the spher- 
ical infall model and the estimation of the collapsed time 
scale have been mainly discussed. Another author evaluates 
the weakly non-linear correction of the background density 
field using the conditional probability distribution function 
around the overdensity and studies the influence of the non- 
linear density inhomogeneities to the dark matter density 
profiles (Lokas 1998). 

In this paper, we investigate the influences of the cos- 
mic background density fleld on the spherical infall model 
using the different approach from the previous works. Fo- 
cusing on the mean expansion around the bound system, 
the Gauss' law implies that the interior of the bound object 
can be regarded as the homogeneous overdensity by taking 
the averaging procedure. Hence, we will treat the averaged 
dynamics of the inhomogeneous gravitational collapse. The 
influence of cosmic density fields on the infalling velocity is 
incorporated into the averaged quantities as the backreacUon 
of the growing inhomogeneities. Importantly, the modified 
dynamics can be non-local. In addition to the mass Af of 
the bound object, the evolution depends on the statistical 
property of the density fluctuations. 

In section |^, we flrst argue how the background den- 
sity fleld modifles the non- linear dynamics (|l|). Then, we de- 
rive the evolution equations for overdensity by averaging the 
inhomogeneous Newtonian gravity. The modified dynamics 
can contain the averaged quantities characterizing the back- 
reaction of the cosmic density inhomogeneity. We explicitly 
calculate these quantities in the Appendix A and obtain 
the simple expressions. The main results in this paper are 
equations (|o|) and (^). Using these expressions, we further 
discuss the infiuence of backreaction effect on the infalling 
velocity in section Final section ^ is devoted to the sum- 
mary and discussions. 



2 AVERAGED DYNAMICS OF SPHERICAL 
INFALL MODEL 

2.1 The backreaction of cosmic density field 

We shall start with the basic equations for Newtonian dust 
fluid: 

dp 



Q^+Vxipv) ^0, 
dv , ^ ^ , ^ ^ 

7^ + (« ■ \/x)V = -Vx<J 
OT 

V\(t> = AirGp, 



(2) 



where t is the proper time of the dust fluid and Vx is the 
derivative with respect to the proper distance X from some 
chosen origin. 



The non-linear dynamics of spherical overdensity (|^) is 
embedded in the above system. Consider the homogeneous 
infalling flow, in which the quantities p, v and (j) are given 
by 



p = p{t) 
V = H{r)X 



<^ = -ttG p(r) \Xf 

where R is the proper radius of the spherical overdensity. 
We then introduce the comoving frame {t, x) defined as 



T = t, X = R{t)x. 

Substitution of into (|^) yields 

p + SHp = 0, 3n + SH + AnGp = 0, 



(3) 



(4) 



where (') denotes the derivative with respect to the time 
t. Because the flrst equation implies the mass conservation 
M = {4n / 3) pR^, the homo geneous system (^) is equivalent 
to the dynamics (^. 

When we take into account the cosmic background den- 
sity field, the dynamics of overdensity affects the evolution 
of background inhomogeneity and the growth of inhomo- 
geneities induces the tidal force, which distorts the homoge- 
neous evolution of the infalling fiows. The situation we now 
investigate is that the background fluctuation is not so large 
and the dynamics of overdensity is dominated by the radial 
motion, that is, the dynamics could be almost the spheri- 
cally symmetric collapse. In this case, we can approximately 
treat the radial infalling flows as the homogeneous system 
and the evolution of the overdensity is solely affected by the 
small perturbation of the background density fields. 

Therefore, the problem is reduced to a classic issue 
of the backreaction. That is, together with the quantities 
(p,TL), the dynamics of the homogeneous overdensity can 
be determined by the backreaction of background inhomo- 
geneities through the non-linear interaction. Phenomenolog- 
ically, the backreaction effect can be mimicked by adding the 
non- vanishing terms characterizing the density inhomogene- 
ity in the right hand side of equations (^. Here, we consider 
the self-consistent treatment based on the equations We 
will derive the dynamics of infalling flows including the back- 
reaction effect. Since the evolution of overdensity is affected 
by the spatial randomness of the cosmic inhomogeneity, it 
should be better to investigate the mean dynamics of the 
overdensity, which can be derived from the system (^) tak- 
ing the spatial averaging. 

Let us divide the quantities p, v and 4> into the homo- 
geneous part and the fluctuating part with the zero mean: 

p = p(r) + 5p(r,X), 

V = H{r)X + u{r,X), (5) 
<!> = ^nG p{T)\Xf + 'i>{r,X). 

Substituting into (^), the equation of continuity becomes 

/dSp 

V dt ' ' R 



{p + mp) + msp + ^v^u 



= ~j^V.{5pu). (6) 
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For the Euler equation, taking the spatial divergence yields 
(37i + 3H' + 47rGp) + ^V^(^^+m+^V^^^ 

= -;^V. [(«.V.)n] . (7) 

In equations and the non- linear interaction of the 
background inhomogeneities are expressed in the right hand 
side of these equations. 

The system and (^) can be divided into the homo- 
geneous part and the fluctuating part by taking the spatial 
average over the comoving volume Vr — (47r/3)ro, where 
the radius ro denotes the comoving size of the characteristic 
overdensity. The spatial averaging is defined as 



1 



(8) 



where {■ ■ ■) denotes the ensemble average taken over the ran- 
dom fluctuations Sp and u. Then, we obtain 

1 



p + 3Hp ■ 



R 



((V. [Spn])), 



R inG 



R 3 37?2 
for the homogeneous part and 

^+SH5p + ^V^u = 0, 
ut H 



(9) 
(10) 



du 

St ' R 
= AnGR^&p, 



(11) 



for the fluctuating part with the zero mean. 

The results ^ , ( [lo| ) and are the desirable forms for 
our purpose. The mean dynamics of overdensity is affected 
by the the averaged quantities {(■ • ■)) corresponding to the 
backreaction of the cosmic density inhomogeneities. These 
quantities are evaluated self-consistently by solving the evo- 
lution equations (|l^ under the knowledge of homogeneous 
quantities p and R. 

2.2 The averaged quantities 

We proceed to analyze the mean dynamics ^ and (|lo|). 
In general relativity, the quantitative analysis of the back- 
reaction terms is still hampered by a more delicate issue 
of the gauge conditions even after deriving the basic equa- 
tions (Buchert 1999). However, in the Newton gravity, we 
are not worried about the evaluation of the averaged quan- 
tities. Let us separate the variables 5p, u and $ into the the 
time- dependent part and the spatial random part: 

Sp = D{t) 5{x), 

u = G R{t)V{t) V:,K{x), (12) 
$ = G R^{t)D{t) k{x), 

where the variables with {") denote the random field. The 
ensemble average is taken with respect to these variables. 
Then, the evolution equations (^, ( p^ and ( pT| ) become 
the four set of ordinary differential equation: 



p + mp 

R i-KG _ 
R + 



= -G(W) Ci, 



(13) 
(14) 



D + 3nD 
V + 2HV 

where 



4nG p V, 
D, 



(15) 
(16) 



Ci 

C2 



((V:. [Spn])) 
G (RVD) 

({Vo; [(«■ Vo,; 

G^iRVy 



are merely the numerical constants. Provided the statisti- 
cal feature of the spatial inhomogeneities S and A, we can 
evaluate Ci and C2. Notice that both 5 and A are related 
through the Poisson equation: 



A: 



4tt S. 



(17) 



Thus, the non-local tidal effect induced by the cosmic den- 
sity fields can be treated in our prescription. 

To evaluate the averaged quantities Ci, C2, it is suffi- 
cient to know the second order statistics of the density in- 
homogeneities, i.e, the power spectrum P{k). Usually, P{k) 
is given by 



{S{k)S{k')) = {2nf5D{k + k') P{k) 



(18) 



in the Fourier representation. However, the naive compu- 
tation using the above definition gives the monopole term, 
i.e, 1 = mode of the spherical harmonic expansion for the 
fluctuations. In our approach, the coherent radial motion 
arising from the monopole contribution has already been 
renormalized in the homogeneous system, that is, the fluc- 
tuations deflned in (jsj) have only the higher multipole 1>1. 
This means that the alternative definition subtracting the 
monopole term 

{5(k)5{k')) = (27r)' P(k) [fo(fc + fc') - — W<5i5(fc - fc')1 (19) 

should be used in our formalism. 

The calculation using (jl^) is essential to obtain the 
non-vanishing averaged quantities. In the Appendix A, the 
details of calculation for averaged quantities are described. 
Here, we present the final expressions: 



Ci 



327r 



6('-o)}. 

{-Co(ro)+3ei(ro)-2C2(ro)}, 



(20) 
(21) 



The function £,i{r\)) denotes the statistical quantity given by 



C;(ro) 



dr 



dk (krf P(k) [ji(kr) 



(22) 



where ji{x) is the spherical Bessel function for l-th order. 

Equations ( |2o| ) and ( |2l| ) are the main results of this 
paper. From these expressions, we immediately see that the 
backreaction effect becomes negligible on large scales due 
to the suppression factor 1/14 . This is correct in the cases 
with the power-law spectrum P{k) oc fc" within the range 
— 3 < n < 1. In the limit ro — > 00, the dynamics ^ and 
( lOl) can be well-approximated by the spherical infall model 
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Figure 1. The averaging quantities Ci,2 as the function of n 
for a power-law spectrum P{k) oc fc". The thin and thick lines 
respectively denote the quantities Ci and Co . which are evaluated 
from Q and using the expressions ( p^ and (ps]). 



3 EFFECTS ON THE INFALLING VELOCITY 

Provided the power spectrum P{k) of the background den- 
sity fluctuations, we can investigate the mean infalling flows 
around the density peak. For the present purpose, it is useful 
to consider the scale-free spectrum in more detail: 

P{k)=Ak". (23) 

Then the quantity ^;(ro) can be expressed in terms of the 
gamma functions. We have 

2"tv r{-n-l)r{l + {3 + n)/2) 



^liro) = ro" — 



n [r(-n/2)]2r(Z-h(l-n)/2)' 



(24) 



The power spectrum is related to the variance of density 
fluctuations ctth, where the subscript th means the vari- 
able with the top-hat smoothing. Identifying the smoothing 
radius with the comoving radius of overdensity ro , the influ- 
ence of the background inhomogeneities to the overdensity 
can be quantified by (Tth. Then, the normalization factor A 
is given by 



47r 



1 2 
(Jth 



[r((2-n)/2)]2r((3-n)/2) 



(25) 



2"-ir(l -n)r((3-Hn)/2) ' 

Thus, the averaged quantities can be characterized by the 
three parameters, ro, n and ctth. 

In Fig.l, we plot the averaged quantities as the function 
of the spectral index n. The typical parameters are respec- 
tively chosen as follows: (Tth = 0.2 and ro = 20 h''^ Mpc 
(solid line); ctth = 0.1 and ro — 30 Mpc (long- dashed 
line) ; CTth = 0.05 and ro = 50 h"^ Mpc (short-dashed line). 
Fig.l says that as the variance (Jth decreases and/or the 
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Figure 2. The ratio of the infalling velocity to the Hubble flow as 
the function of the spectral index n (long-dashed, short-dashed, 
dotted, dot-dashed-lines). By solving the mean dynamics (^^- 
( [l6| ) , the quantity given by (^) is evaluated at the time a/ /oi = 4 
in the Einstein-de Sitter universe. We choose the initial density 
contrast <5 = 0.3. The initial condition for the fluctuating part 
are described by the linear growing mode of perturbations ne- 
glecting the backreaction effect. The results are compared to the 
homogeneous infalling flow obtained from the simple dynamics 
(|l|) (4/im-so/jd line). 



radius ro increases, the quantities Ci and C2 become neg- 
ligible. The thick lines shows the variable C2. Interestingly, 
C2 turns its signature. The negative value is obtained at 
n < — 1, while we have the positive value C2 at n > — 1. By 
contrast, the quantity Ci denoted as the thin lines becomes 
positive and almost constant over the range —2.8 < n < 0.8. 

Fig.l indicates that backreaction of the cosmic den- 
sity inhomogeneities can prevent the spherical shell from 
infalling into the center of overdensity for the negative spec- 
tral index (see eqs.[^[Q). The effect becomes noticeable 
when n < —1. The result can be interpreted as follows. 
The quantities Cais expressed in terms of the geometri- 
cal optics. From (|l2), we can write the peculiar velocity as 
diUj = (6/3)5ij -f aij, where 6 and aij are the expansion 
scalar and the shear tensor, respectively. Then we have 

({V. [(u ■ V^)u] )) = ((ie^)) + + {{u ■ Ve». (26) 

The last term in the right-hand side of ( |26| ) means the flows 
of inertial frame. It is easy to show that the first and the 
second terms in the right-hand side of ( p^ always become 
positive. On the other hand, the last term can be written by 



G2(J?V)2 



327r 



^(2/ + l)Ci(ro)-ei(ro) 



(27) 
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which is always negative. This means that the anisotropic 
flows of cosmic density field can induce the effective pres- 
sure and this could dominate the expansion and the shear 
terms. The effect becomes significant when the long wave- 
length fluctuations have the large power, corresponding to 
the small spectral index n < —1. Accordingly, the back- 
ground inhomogeneity apparently prevents the collapse of 
the averaged overdensity. 

To see the backreaction effect more explicitly, we ana- 
lyze the infalling velocity. The flow field around the spherical 
overdensity can be quantified by 



H\r\ 



H ~n 

H ' 



(28) 



The above quantity can be evaluated by solving the evolu- 
tion equations (p^)-([l6|). 

In Fig. 2, we plot the ratio of the infalling velocity to the 
Hubble flow as the function of the spectral index. Here, the 
cosmic expansion is assumed to be that of the Einstein-de 
Sitter universe, a = (tjtif'^^ . The initial velocity flow 7i of 
the overdensity is obtained by equating it with the Hubble 
expansion. As for the initial overdensity, p can be character- 
ized by the density contrast & — p/ pt — 1, where pb is the 
homogeneous density of the universe. We choose the initial 
density contrast as5 — 0.3. On the other hand, the evolution 
of the fluctuating part D and V can be initially described 
by the linear growing mode neglecting the backreaction ef- 
fect. Then, we solve the evolution equations by varying the 
spectral index n. The results are depicted when the final ex- 
pansion factor a/ becomes four times larger than the initial 
expansion factor m, i.e, a/ /ai = 4, within the validity of the 
approximation neglecting the higher order perturbations. 

In Fig. 2, we also plot the ratio Vini/H\f\ evaluated from 
the spherical infall model (|l]) {thin-solid line). Clearly, the 
deviation from the spherical infall model becomes signifi- 
cant for the small radius of the overdensity and the large 
variance of the fluctuations {thick-long- dashed, thick-short- 
dashed lines). The cosmic background inhomogeneity can 
weaken the infalling velocity. This effect leads to the over- 
estimation of the cosmic expansion when we compare the 
observation of the infalling velocity with the spherical infall 
model. Consequently, we would prefer the low density uni- 
verse. However, the deviation becomes small as the spectral 
index increases. The backreaction effect becomes negligible 
and the infalling velocity can be well-described by the spher- 
ical infall model. At the time a; jai — 4, the index n ~ 0.5 
is marginal, which is rather larger value than n — —1. 

As another example, a realistic case with the cold dark 
matter (CDM) spectrum is examined. In Fig.3, we specifi- 
cally consider the standard CDM, i.e, f2o = 1 and Ao = 0. 
The fitting form of the power spectrum given by Bardeen et 
al. (1986) is used and normalized by the top-hat fluctuation 
amplitude at 8/i~^Mpc, erg = 0.6, according to the cluster 
abundance (Kitayama & Suto 1997). The Hubble param- 
eter is given by Ho = 100 /ikms~^/Mpc with h = 0.5. In 
Fig.3, the present values of Ci,2 (upper-panel) and Vinf/fflrl 
(lower-panel) are depicted as the function of the averag- 
ing scale ro. The initial conditions for the lower-panel are 
set at 2 = 5, similar to the Fig. 2. The figure states that 
both quantities Ci and C2 become positive on large scales 
ro > 30/i~^Mpc and the resulting ratio Vinf/-ff|r| takes the 
slightly lower value, compared to the original spherical in- 
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Figure 3. Scale-dependence of the averaging quantities Ci_2 
(upper-panel) and the ratio of the infalling velocity to the Hub- 
ble flow (lower panel) in the case of the standard CDM model. In 
both panels, we used the power spectrum P{k) given by Bardeen 
et al. (1986) normalizing at the present epoch with the parame- 
ters (f2o, Ao, /i, eg) = (1.0,0.0,0.5,0.6). Upper-panel: the averag- 
ing quantities Ci {thin-solid line) and C2 {thick-solid line) eval- 
uated at present ; Lower-panel: the present values of the ratio 
Vinf/^^M evolved from 2 = 5 for the original spherical infall 
model {thin-solid line) and our model {open triangles). 



fall model {thin-solid line in lower-panel). This is because the 
dominant contribution of the CDM spectrum to the averag- 
ing quantities is the low-fc part with the effective spectral 
index — 1 < UeB < 0, where rieff = dlog P(fc)/dlog (see 
Fig.l). The behavior of infalling velocity is thus understood 
from Fig. 2. On the other hand, the deviation of Vinf/iif|r| 
from the spherical infall model becomes large as decreasing 
the averaging size ro. The quantity C'2 turns its signature 
at the relatively small size, ro ~ 9/i~^Mpc. In any cases, 
the present model indicates that the method using infalling 
velocity still leads to the underestimation of the density pa- 
rameter in real universe. 



4 DISCUSSION AND CONCLUSION 

We have analyzed the mean dynamics of the spherical infall 
by taking into account the effects of cosmic background den- 
sity field self-consistently. After deriving the evolution equa- 
tions, we have rigorously computed the averaged quantities 
appearing in the equation of the homogeneous overdensity. 
The main result in our analysis is equations ( po| ) and (^l[). 
The expressions indicate that the long-wavelength inhomo- 
geneity with the larger power might induce the anisotropic 
flows and this could dominate the expansion scalar and the 
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shear tensor in the case of the negative spectral index for 
the scale-free power spectrum. Thus, the infalling velocity 
around the averaged overdensity could be evaluated as the 
rather small value, compared to the spherical infall model. 
Solving the evolution equations we have confirmed 

these things. For n < 0, the infalling velocity significantly 
deviates from the prediction of the spherical infall model. 
The same result is obtained in the case of the CDM spec- 
trum for which the dominant contribution of the background 
inhomogeneity comes from the spectrum with the effective 
index riefi < 0. This indicates that we underestimate the 
density parameter Q when using the simple spherical infall 
model. It might lead to another suggestion that the esti- 
mation of the cluster abundance using the Press-Schechter 
formalism is changed when we take into account the effect 
of the cosmological background density field. On the other 
hand, for the spectral index n > 0, the effect of background 
inhomogeneity is negligible and the simple spherical infall 
model could be a good approximation for the infalling ve- 
locity. 

The analysis of the infalling velocity with n < is quali- 
tatively consistent with the numerical simulation of the clus- 
ters of galaxies. Lee et al. (1986) analyzed the N-body model 
of a flat universe with an initial Poisson distribution and 
concluded that the prediction of the spherical infall model 
underestimates the density parameter in the mean though 
the scatter is large. Hoffman (1986) presented an analytical 
model of a spherical collapse with a global shear. He showed 
that when the ratio of the anisotropic flow to the isotropic 
infall velocity exceeds about 50%, infall velocities are higher 
than in the pure spherical infall model. In other cases with 
the small anisotropy, corresponding to the situation consid- 
ered here, the simple dynamics (^) underpredicts the density 
parameter, which is consistent with our result. Villumsen & 
Davis (1986) studied the nature of the velocity held around 
large clusters in SI = 1 cosmological N-body simulations 
with more realistic initial conditions. They concluded that 
the substantial subclustering seen on small scales causes a 
systematic underestimate of the density parameter. Lilje & 
Lahav (1991) also considered the averaged density and infall 
velocity profiles around clusters of galaxies in biased CDM 
models and showed that the prediction of the spherical in- 
fall model overestimates the density parameter in the case 
of the bias parameter 2.1. As authors themselves stated, the 
several reasons could cause the discrepancy: the selection 
and the definition of the clusters, in connection with the 
normalization and the bias parameter in the N-body sim- 
ulation, and/or another interpretation that the large-scale 
shear will make the infall velocity larger. It is therefore dif- 
ficult to compare their result with ours. Here we adopt the 
former reason which would be most likely and conservative. 
Then, their result is reconciled with ours provided that the 
bias parameter is taken to be less than 1.8. Thus, all of the 
other works could be in agreement with the present analysis 
and provide us a simple physical picture: when the spectral 
index n < 0, the anisotropy of the random velocities can 
act as an effective pressure and yield the slowing down the 
infalling velocity. 

Although we could not flnd any simulations in the cases 
with n > 0, the similar behavior can be found from the evo- 
lution of the power spectrum (Makino, Sasaki & Suto 1992, 
Jain & Bertchinger 1994). The weakly non-linear analysis 



using the perturbation theory shows that the non-linear evo- 
lution of the power spectrum has the different growth rate, 
depending on the spectral index n. For the single power-law 
spectrum n < — 1, the non- linear growth of power spectrum 
significantly enhances, while the growth is suppressed if we 
have the index n > — 1. These behaviors remarkably coincide 
with our calculations of the averaged quantities. Therefore, 
we might expect that our result remains correct even in the 
n > case. The validity of our prediction will be confirmed 
by the numerical simulations elsewhere. 

In this paper, we have found the significant contribu- 
tion of the background inhomogeneities to the infalling flow 
around the overdensity. However, the influences of cosmic 
background density field on the local bound system deserve 
further study. The background inhomogeneities also affects 
the internal structure of the bound object. The formation 
of substructure near the central region significantly alters 
the dynamical evolution. The non-radial motion induced by 
the tidal torques might slow the dynamics of the collapse, 
which leads to the previliarization process (Peebles 1990). 
We think that the background density inhomogeneity could 
be one of the most important sources of the previliarization 
effect. To resolve this issue, the further analytical study must 
be developed. 
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APPENDIX A: CALCULATION OF THE 
AVERAGED QUANTITIES 

In this appendix, using the definitions (|) and (|l|), we cal- 
culate the averaged quantities appearing in the right hand 
side of ([o|) and (|). 

Using the new variables (p^, the averaged quantities 
can be written as 



Ci = A, 

Ca = -(47r) A-{^nfB^C, 



(Al) 
(A2) 



where we used the relation (^) . The quantities A, B, and C 
are respectively given by 

A = ({a, [s d^A] )), B = iP)), c = ({{d.d.kwd'k))), 

where di and 9' are the derivative with respect to and 
the summation convention is used. In terms of the Fourier 
representation, the above equations become 



A 



47r / 3 



(27r)6 



1 



2{kik2) 



<fki(fk2 

(27r)6 



X {5{ki)5{k2)) e-'<'=i+'=^'"^ 
{5{ki)5{k2)) e-'(*^i+*^2'", 



(A3) 
(A4) 



B = ^7 V(2; + l)e,(ro) 

TV Vr ^-^ 



c = 



327r 

IT 



oo 

^(2; + l)C;(n.)-ico(n,)-^C2(ro) 



3^^ ' 3 

where the quantity ^; is defined by (22p. Substituting the 
above expressions into and (|A2), we finally get ( |20[ ) 

and 



(rkid'^k2 f kl ■ k2 



(A5) 



We first perform the volume integral Vr using the following 
formula : 



d X e 



oo L 
1=0 m = — l 

r° 2 

X / dr r ji{kir)ji{k2r), 
Jo 



where Yim is the spherical harmonics and the argument fli 
denotes the angular position for the wave vector ki. 

To proceed further, we rewrite the integral ki with the 



spherical coordinate s (k j , 9i 



Then, we substitute the 



definition ( |19[ ) into (A3), (A4) and (AS) and evaluate the 
integral over 9i and (j>i- For A and C, the calculation is 
tractable if we use the relations : 



kl ■ k2 



\ k\k2 J 
where we define 



47r 



Fo(!^l) • Yo(Sl2) + |?2(!^l) ■ Y2(Sl2) 





Yl{Q.l)-Yl{Q.2)= YCm{^l)Ylm{0.2 



(A6) 



With the careful manipulation, the following results can be 
obtained: 

A = l.{i,,[ro)-Uro)}, 
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